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ABSTRACT 

We present the initial results from the Spitzer Survey of the Small Magellanic Cloud (S 3 MC), 
which imaged the star-forming body of the Small Magellanic Cloud (SMC) in all seven MIPS and 
IRAC wavebands. We hnd that the Fg/F24 ratio (an estimate of PAH abundance) has large spatial 
variations and takes a wide range of values that are unrelated to metallicity but anticorrelated with 24 
/xm brightness and F24/F70 ratio. This suggests that photodestruction is primarily responsible for the 
low abundance of PAHs observed in star-forming low- metallicity galaxies. We use the S 3 MC images 
to compile a photometric catalog of ~ 400,000 mid- and far-infrared point sources in the SMC. The 
sources detected at the longest wavelengths fall into four main categories: 1) bright 5.8 /im sources with 
very faint optical counterparts and very red mid-infrared colors ([5.8] — [8.0] > 1.2), which we identify 
as YSOs. 2) Bright mid-infrared sources with mildly red colors (0.16 < [5.8] — [8.0] < 0.6), identified 
as carbon stars. 3) Bright mid-infrared sources with neutral colors and bright optical counterparts, 
corresponding to oxygen-rich evolved stars. And, 4) unreddened early B stars (B3 to 09) with a large 
24 fjm excess. This excess is reminiscent of debris disks, and is detected in only a small fraction of 
these stars (< 5%). The majority of the brightest infrared point sources in the SMC fall into groups 
one to three. We use this photometric information to produce a catalog of 282 bright YSOs in the 
SMC with a very low level of contamination (~ 7%). 

Subject headings: dust — galaxies: dwarf — galaxies: stellar content — infrared: ISM — infrared: 
stars — Magellanic Clouds 



1. INTRODUCTION 

The Small Magellanic Cloud (SMC), one of the clos- 
est and most prominent neighbors of the Milky Way, 
is a southern hemispher e dwarf galaxy of low mass 
(M dyn ~ 2.4 x 10 9 M Q ; iStanimirovic et all 120041) and 
small size (R* ~ 3 kpc). The SMC is actively interacting 
with its companion the Large Magellanic Cloud (LMC, 
studi ed by the Spitzer SAGE project; iMeixner et al.l 
2006). Both galaxies are being tidally disrupted by the 
Milky Way, and also suffer from ram-pressure interac- 
tions with the halo of hot gas s urrounding the Galaxy 
(e.g.. iGardiner fc Noguchil 119961) . The results of this 
complex dynamical interplay are the dramatic features 
of the Magellanic Bridge and Magellanic Stream, which 
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dominate the high-velocity H I emission in the Southern 
sky. 

Befitting its nature as a gas-rich late-type dwarf galaxy, 
the SMC contains a number of sites of active star forma- 
tion that are apparent in Ho images, although none are 
as prominent as the 30 Doradus star burst in the LMC. 
The brightest SMC star forming complex is the giant 
H II region N 66 (also known as DEM 103, or NGC 346; 
lHenizeill956t iDavies e t al.lll97 6j). which contains approx- 
imately 60 O stars l[Massev et al.lll989l) . The star forma- 
tion activity of the SMC, as well as most of its gas, are 
concentrated in the feature known as the "Bar" (proba- 
bly the body of the SMC proper). A second feature, the 
"Wing" , is an extension roughly perpendicular to the 
Bar on its eastern side that continues into the H I bridge 
that joins both Clouds (Figure^. As one may expect 
in an interacting system, the three-dimensional struc- 
ture of the galaxy remains quite unclear. Recent anal- 
yses show that the prominent irregular features, which 
are apparent in the light of young star s, disappear when 
looking at the old stellar populations |Cioni et aIll200Ct 
iZaritskv et ail 120001 iMaragoudaki et all I200i|k These 
are regularly distributed, and the shape of the SMC 
appears spheroidal. Observations indicate that the dis- 
tances to the Bar and the Wing are measurably differ- 
ent, and that the line-of-sight depth of the SMC is not 
negligible, amounting to 6 to 12 kpc <)Crowl et al.ll200ll 
and references therein). The distance to the center of 
the SMC is approximately 61.1 kpc, corresponding to 
a distance modulus of m — M = 18.93 mag and a spa- 
tial s cale of 0.3 pc/" i|Hilditch et al.l2005llKeller fc Woodl 
2006). To put this in context, the spatial resolution 
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of the 3.6 — 8.0 /an data presented here is similar to 
the 60 /tm resolution of the Infrared Astronomical Satel- 
lite (IRAS) at 2 kpc. The thickness of the SMC intro- 
duces a dispersion of ~ 0.2 mag in the distance mod- 
ulus for objects located near either edge. The line-of- 
sight to the SMC is very clear, with a foreground ex- 
tinction Ay ~ 0. 12, corresponding t o a reddening of 
E{B - V) ~ 0.04 llSchlegel et al.lll99l . 

The proximity of the SMC makes it the lowest metal- 
licity gas-rich galaxy that can be studied in great de- 
tail. For observational reasons metallicity, Z, is com- 
monly quantified using the oxygen abundance relative to 
hydrogen. The gas-phase oxygen abundances measured 
in SMC H II regions are 12 + log(0/H) » 8.0 lIDufourl 
Il975l IDufour k HarlowllT977l IDufourl 1 1984ft . only a fac- 
tor of 4 to 5 higher than those observed in the lowest 
metallicity galaxies known (I Zw 18 and SBS 0335-052), 
and a factor of 5 lower than solar metallicity. The same 
metallicity with very little dispersion is observed in sev- 
eral stellar types in the SMC, in cluding K, F, and A su- 
pergiants and Cepheid varia bles llRussell &: Bessel]Hl989t 
IHmiil997t ILuck et al.lll998t IVenrJ 1 1999ft . Its low heavv- 
element abundance suggests that the SMC may be one of 
the best local templates for studying the primitive inter- 
stellar medium and providing insight into the processes 
at work in primordial galaxies at high redshift. In par- 
ticular, the interaction between massive stars, gas, and 
dust is unusually accessible in this nearby galaxy. 

Compared to the Milky Way, the dust properties of 
the SMC are peculiar. For quite a while it has been rec- 
ognized that the UV ex tinction curve in th e SMC lacks 
the 2175 A"hump" (e.g.. lGordon et al.l2003(l which is in- 
terpreted as signaling the lack of a populatio n of small 
carbonaceous dust grains l|Li fc Drainel ^002). Indeed, 
most current models of dust composition in the SMC 
rely heavily on silicates to reproduce its extinction curve, 
with typical silicate to carbon ratios by m ass in the range 
2-12 fe.g.. lWeingartner & Dramell2001j) . Whether this 
dearth of carbonaceous grains is maintained for the large 
grains remains unclear, as there is a well known anticor- 
relation between the specific frequency of carbon stars 
(the main producer s of carbonaceous g rains) and galac- 
tic metallicity (e.g.. lBlanco et aLlll980ft . that results in a 
large abundance of these stars in the SMC. 

Other lines of research also suggest a dearth of the 
smallest dust grains in the low-metallicity interstel- 
lar medium of the SMC. The first wide-field infrared 
maps of t he SMC were acqu i red b y I RAS and pre- 
sented by iSchwering &: Israel l)1989f) . iSauvage et alJ 
(1990) studied the IRAS colors of the SMC, compar- 
ing them to the emission from larger, higher metallic- 
ity galaxies, and concluded that there was evidence for 
a paucity of the small grains that are responsible for the 
bulk of the extended 12 /im emission. More recently, 
IStanimirovic et all l)2000ft reprocessed the IRAS data to 
produce maps with improved angular resolution of 1' at 
12 /im, 25 /zm, and 60 /an, and 1.'7 at 100 /an. By com- 
paring these data with H I observations of the SMC at 
similar resolution, the y were able to study the large-scale 
properties of the ISM. IStanimirovic et all used these far- 
infrared (FIR) observations to measure the dust mass, 
the dust temperature and its spatial variations, and the 
dust-to-gas ratio, DGR, across the SMC. This study 
found a mean dust-to-gas ratio in the SMC that is a 



factor of ~ 100 times lower than in the Milky Way, pro- 
viding a different environment for star formation to take 
place. It also concluded that the dust in the SMC is dom- 
inated by large grains, with weaker contributions from 
polycyclic aromatic hydrocarbons (PAHs, typical sizes 
4A< a <12A; Des ert et a.l.lllQQO^ and very small grains 
(VSGs, 12A< a <150A) than are found in the Galaxy. 
One of the limitations of this study is the use of the 
60 /im IRAS waveband to determine dust temperature. 
This waveband is expected to have a substantial contri- 
bution from stochastically heated small dust grains, out 
o f equilibrium with the radiation field. 

iBot et alJ 1)20041) expanded this work by adding In- 
frared Space Observatory (ISO) 170 /im imaging of the 
SMC to the IRAS data. Using a combination of 100 and 
170 /im data they measured an appreciably higher dust- 
to-gas ratio, ~ 30 times lower than Galactic. They also 
modeled the spectral energy distribution to decompose 
the emission into the contributions by PAHs, VSGs, and 
large grains, finding a substantial 60 /tm excess. Such 
excess could be caused by an enhanced interstellar radi- 
ation field in the SMC, or perhaps by a change in the 
grain size distribution with respect t o the Galax y, with 
more VSGs relative to large grains. IBot et all suggest 
that grain destruction, perhaps shattering by supcrnovae 
shocks, is the driver behind the different properties of the 
dust in the SMC and the Milky Way. 

These studies point to a dust-to-gas ratio in the SMC 
much lower than Galactic, scaling with its metallicity as 
DGR oc Z 2 , a res ult similar to that found f or a collec- 
tion of galaxies bv lLisenfeld &: Ferraral l)1998ft . Note that 
measurements of the mass-loss rate of the most evolved 
asymptotic giant branch (AGB) stars in the SMC, how- 
ever, find a much larger du st-to-gas ratio, scaling as 
DGR oc Z 1 l|van Loonl 12006ft . Because these stars are 
believed to dominate the ISM dust production, this ob- 
servation suggests that efficient destruction of dust takes 
place in the metal-poor star-forming ISM, along the lines 
of the evidence mentioned in the previous paragraphs. 

Detailed modeling o f SMC UV absorpti o n an d 
FIR emission data by iWeingartner fc Drainel (2001), 
iLi fc Drainel {2002), and iClavton et all 1)20031) indicate 
that a deficit of PAHs, as well as the lower abundance of 
VSGs suggested by IRAS, and perhaps a change in com- 
position with more silicate and fewer gra phite grains are 
requir ed to reproduce the observations. iContursi et alJ 
(2000) studied N 66 using the cameras onboard ISO, 
and found that the emission from the Aromatic Infrared 
Band (AIB) in the 6.75 /im ISO waveband, commonly 
attributed to PAHs, was faint compared to the 15 /im 
dust continuum. This was confirmed spectroscopically 
using ISO's continuously variable filter (CVF) to ob- 
tain low-resolution spectra toward emission peaks, all of 
which show v ery weak or nonexis tent AIB emission. As a 
counterpoint, iReach et all l)2000D used the CVF to obtain 
spectroscopy toward a quiescent molecular cloud in the 
SMC. They measured AIB-to-continuum ratios similar 
to those observed in the Milky Way, pointing to similar 
relative abundances of PAHs. These data suggest that 
PAHs do form and survive in the SMC and other star- 
forming low-metallicity galaxies in quantities very simi- 
lar to Milky Way PAHs, where the environment allows it. 
Indeed, SMC observations show us that dust properties 
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are environment dependent, and one of the main moti- 
vations for a new infrared survey of the SMC is to study 
and characterize the influence of the local environment 
on the properties and composition of interstellar dust. 

The Spitzer Space Telescope allows dramatic improve- 
ments in our ability to study the ISM of the SMC, with 
previously unattainable images at 3.6 to 8.0 /zm and a 
factor of 4 — 10 higher angular resolution and typically 
an order of magnitude better sensitivity at FIR wave- 
lengths than previous missions. The Spitzer data enable 
us to take a giant step forward in the study of the de- 
tailed, distribution of all types of dust grains, including 
PAHs which primarily emit at wavelengths shorter than 
12 /im. The high angular resolution of Spitzer means 
that point sources can be separated from diffuse emis- 
sion out to 70 /im, allowing investigations of large sam- 
ples of IR point sources (such as carbon stars and young 
stellar objects) and simultaneously providing a cleaner 
picture of the ISM. The unprecedented sensitivity and 
resolution of the Spitzer data makes it possible to deter- 
mine the abundance and distribution of dust grains and 
PAHs in relationship to their environment, to test and 
constrain dust models, and to obtain a complete census 
of star formation in the SMC, the nearest "primitive" 
galaxy. 

In this paper we present an overview of the Spitzer Sur- 
vey of the Small Magellanic Cloud (S 3 MC), discussing 
the observations, methodology, and some of the initial 
results on the extended emission and the infrared point 
sources in the SMC. This is the second in a series of 
papers presenting the r esults of the survey (see also 
iStanimirovic et all 2005). In |JS|we describe the obser- 
vations, in Sj3 we describe in detail the data reduction 
process employed to obtain the final S 3 MC images, in 
21 we discuss the point source photometry for the global 
survey and the comparison with existing datasets, in iJ3 
we analyze the extended dust emission and discuss the 
different stellar populations detected by Spitzer, and in 
SjJJwc summarize our findings. 

2. OBSERVATIONS 

We obtained Spitzer Space Telescope images of the 
Small Magellanic Clou d in all the bands of the Infrared 
Array Camera (IRAC. lFazio et al.ll200l and the Multi- 
band Imaging P hotometer and Spectrometer (MIPS, 
iRieke et al.l 12004]) as part of the Cycle I General Ob- 
server project 3316, "The Small Magellanic Cloud: A 
Template for the Primitive Interstellar Medium". This 
project was designed to map as large an area of the SMC 
as possible with good sensitivity and within the scope 
of a small Spitzer proposal (less than 50 hours), with 
the goal of sampling a diverse set of environments. The 
images cover an area of approximately 2.8 to 4 deg 2 in 
the various bands, with a region of ~ 2.5° x 1° mapped 
along the SMC Bar, and ~ 1.5° x 1° along the Wing, and 
the Astronomical Observing Requests (AORs) were de- 
signed to provide contiguous coverage without imposing 
scheduling constraints (Figure HJ. 

The MIPS data were acquired in November 2004, using 
five large AORs each lasting over 3 hr, and two smaller 
AORs slightly over 1 hour long each. The angular resolu- 
tion of Spitzer at 24, 70, and 160 /im is 6", 18", and 40" 
respectively. The observations used the MIPS Scan As- 
tronomical Observing Template (AOT), with a medium 



(675 s _1 ) scan rate as a compromise between coverage 
and sensitivity. A total of over 11,700 individual images 
were obtained in each of the three MIPS bands. The 1-a 
sensitivities measured in the MIPS images are approxi- 
mately 0.045 MJy sr" 1 at 24 /im, 0.6 MJy sr" 1 at 70 /mi, 
and 0.7 MJy sr -1 at 160 /im. The sensitivity is degraded 
at 70 /im by the pattern noise present in the image. At 
160 /im it is questionable whether it is possible to find 
regions of the image that are free of emission in which 
to measure the noise level, so the sensitivity estimate 
for this wavelength has a significant uncertainty. Fur- 
thermore, to accommodate the observations within the 
scope of the proposal the oversampling factor (i.e., the 
number of times a given position in the sky is observed 
with different parts of the detector) at 160 /tm was two 
rather than the recommended four. The corresponding 
1-a point source sensitivities are 0.043, 5, and 30 mJy 
respectively. 

The IRAC data were obtained during 2005 May us- 
ing 8 large 4 hour-long and 2 small 2 hour-long AORs 
with a frame time of 12 seconds and a 3 point cycling 
dithering pattern. The point response function of IRAC 
has a FWHM of approximately 1766, 1772, 1788, and 
1798 in the 3.6, 4.5, 5.8, and 8.0 /im bands. The typical 
oversampling factor of a given position in the sky was 
~ 12. Anticipating a number of bright point sources, the 
AOTs for the fields along the SMC Bar were obtained 
using the high dynamic range mode, which takes a short 
0.6 second exposure along with each 12 second exposure 
on the sky. This allows us to recover the correct fluxes 
for very bright stars that are saturated in the 12 second 
exposures. The saturation limits for stars in an 0.6 sec- 
ond exposure are 630, 650, 4600, and 2500 mJy in the 
3.6, 4.5, 5.8, and 8.0 /im bands respectively. Approx- 
imately 6,700 individual images were obtained at each 
band. The 1-a sensitivities estimated from the mosaics 
are 0.008, 0.01, 0.04, and 0.04 MJy sr" 1 , corresponding 
to point source sensitivities of 0.7, 1.1, 8, and 6.5 /iJy 
respectively. Table ^ summarizes the survey properties. 

3. DATA REDUCTION 

We used the MOPEX data reduction package (version 
050905) provided by the Spitzer Science Center (SSC) 
to combine the individual Basic Calibrated Data (BCD) 
frames into large mosaics. There are a number of possi- 
ble methods for processing the data into a final mosaic 
of the SMC. We chose to begin with the BCD frames, 
which are derived from a single data collection event (a 
DCE, or a single frame exposure). The BCD images 
were generated by the SSC using the standard pipeline 
(version S11.4). The detailed pipeline processing of raw 
im ages from Spitz er to BCD FITS files is f ully d iscussed 
bv iGordon eta!] p005]) and iMasci eta!] pDOlj) . Each 
BCD image is accompanied by an uncertainty and a mask 
image, which allow proper weighting of the data and re- 
moval of corrupt pixels. 

The SSC pipeline provides, in addition to the BCD 
images, a preliminary mosaic of each AOR referred to 
as the post-BCD mosaic. We use these images primarily 
for comparison when we generate the full, much larger 
mosaic of the SMC. 

3.1. MIPS Data Reduction 
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Fig. 1. — Spatial coverage of the S 3 MC data, overlaid on the Ha "parking lot camera" image of the SMC (Kcnnicutt et al.lll995l) . 
The white outline indicates the coverage of the IRAC observations, while the gray outline shows the coverage for MIPS. Note that this 
is only approximate, since the placement of the cameras on the Spitzer focal plane results in slightly different areas of the sky mapped 
for the different wavebands. The SMC Bar is the region where the bulk of the star formation (indicated by Ha emission) takes place, 
extending roughly from (01 h 15 m ,-73°30') to (00 h 54 m ,-72°00')- The Wing region, notable in H I but with very little diffuse Ha emission, 
is approximately perpendicular to the bar and extends eastward from (01 h 00 m ,— 72°45') to (00 h 45 m , — 73° 15') and beyond. 



TABLE 1 
Survey Properties 



Property 


Unit 


3.6 £tm 


4.5 (im 


5.8 fim 


8.0 /im 


24 fim 


70 /im 


160 /^m 


Sensitivity 1 


MJy sr _1 


0.008 


0.010 


0.040 


0.040 


0.045 


0.600 


0.700 


Sensitivity 2 


MJy 


0.7 


1.1 


8 


6.5 


43 


5,000 


30,000 


Resolution 


arcsec 


1.66 


1.72 


1.88 


1.98 


6 


18 


40 


Saturation 


mJy 


630 


650 


4,600 


2,500 


1,000 


5,750 


750 


Area 


deg 2 


2.77 


2.77 


2.77 


2.77 


3.94 


3.70 


3.69 


Source Count 3 




292,000 


289,000 


80,000 


62,000 


16,000 


1,800 




Completeness Flux 4 


MJy 


31 


24 


80 


80 


400 


40,000 




Photometric Zero Point 5 


Jy 


277.1 


179.4 


113.9 


63.1 


7.253 


0.8198 





1 For extended emission, 1-cr. 
2 For point sources, 1-cr. 
3 Approximate. 

4 Approximate. See discussion in H4.2I 

5 Martin Cohen, priv. comm.. Consistent with Reach ct al. 12005). 



Inspection of the 24 /im MIPS BCD and post-BCD 
images revealed a small number of saturated sources and 
obvious latent images from bright sources or cosmic ray 
hits. No attempt was made to explicitly remove cos- 
mic rays from individual BCD images before combining 
frames together. To avoid including cosmic ray hits in 
the final mosaics, we used the redundancy and outlier 
rejection in mosaicking. Only temporal outlier rejection 



was performed, where redundant fields taken at different 
times are compared to remove artifacts. 

The 70 /im data contain several instrumental effects as 
well as bright extended emission that make point-source 
extraction difficult. The SSC provides two types of BCD 
images: normally processed (unfiltered) and those with a 
time-median filter applied that removes most of the back- 
ground signal (filtered). Following the recommendation 
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of the SSC, we used the 70 /im filtered data to obtain 
point-source fluxes. However, only the unfiltered data 
preserve extended emission. In this paper, we present 
mosaics made from the unfiltered 70 /im data. Because 
of the complexities of dealing with the "first frame" ef- 
fect, row-droop, a nd other similar artif acts in the current 
software pipeline l|Gordon et al.ll2005j) . we have not ap- 
plied any extensive processing to form the 70 /im mosaic. 

We combined the individual MIPS observations to 
form a mosaic using the default linear interpolation 
method. Adjacent, partially overlapping fields were com- 
bined using corrections computed by the MOPEX task 
overlap.pl, which matches fluxes in common regions to 
guarantee a consistent background level. Because this 
process is computer intensive, we applied it to subsets of 
~3,000 images at a time. MOPEX calculates the offsets 
corresponding to each image by generating a matrix of 
median values for different DCEs and solving a system 
of linear equations. 

We linearly combined the background-corrected MIPS 
BCD images, weighting them by the number of cover- 
ages at a given position and the associated uncertainty 
image. The outlier rejection settings employed are de- 
scribed above. 

3.2. IRAC Data Reduction 

For the IRAC data, we performed a number of pre- 
processing steps on the individual BCDs before com- 
bining them into a complete mosaic. We used the 
full "dual-outlier" rejection method (see following para- 
graphs) in order to address bright, solar system objects 
whose proper motion is easily detectable during our dwell 
time on the SMC. Following recommendations from the 
SSC (D. Makovoz, priv. comm.) we used the default 
linear interpolation method to resample pixels from the 
BCD images into the new fiducial mosaic image frame in 
all instances. 

Following this correction, we used specialized MOPEX 
routines to remove two artifacts present in the IRAC im- 
ages: the so called muxbleed and column-pulldown cos- 
metic effects caused by the presence of bright sources in 
individual DCEs. Muxbleed manifests itself as a trail 
of repeating bright pixels (every fourth pixel) along the 
readout direction of the chip after detection of an ex- 
tremely bright star. The effect of column-pulldown is 
to bring down the co unts in an entire co lumn associated 
with a bright source l|Patten et alJl200"4j) . These artifacts 
are dealt with using the MOPEX script cosmetic.pl. 
Pixels are masked out of the final image processing if 
both 1) at least one pixel is found to exceed 50 MJy/sr 
and 2) if more than 10% of the pixels along the bright 
pixel row are found to have fluxes outside of the range 
2-40 x the RMS value for the row. These settings are 
the default values, found to work best in removing the 
muxbleed problems. 

The next pre-processing step removed the saturated 
pixels in those fields (along the Bar) for which we ob- 
tained high dynamic-range imaging. We compared the 
short and long exposures for each pixel and identified sat- 
urated pixels. Those pixels found to be saturated in the 
long exposure were then replaced with the corresponding 
pixel in the short exposure. 

An important component of the mosaicking of the 
IRAC BCDs is outlier rejection. We combined the cor- 



rected IRAC BCD images by performing both time- 
domain (outlier rejection) and the more elaborate 
spatial-temporal (dual-outlier rejection) filtering pro- 
cesses. We used the recommended (default) settings for 
the values of the relevant modules in outlier rejection. 
We linearly interpolated the BCD data onto the final 
mosaic. The images are combined as a weighted sum, 
using the associated uncertainty images, to form the fi- 
nal image. 

The IRAC images from the 5.8 /im waveband (and, to 
a lesser extent, those at 8.0 /im) presented an additional 
challenge in the form of large-scale additive gradients in 
the background throughout the mosaic. The gradient is 
caused by the combination of a time varying offset in the 
readout amplifier, and the mosaicking algorithm imple- 
mented in MOPEX (SSC, priv. comm.). Following the 
recommendations of the SSC, we made a "super dark" 
frame using the BCDs within each AOR. Using proce- 
dures developed in IDL, we stacked the BCDs for each 
AOR and computed a pixel-by-pixel median value. We 
then subtracted the overall median of this super dark 
from each pixel, so that its DC level was zero. We sub- 
tracted the super dark from each BCD, flat fielded the 
images, and continued with the same mosaicking routines 
that were used for the 3.6 and 4.5 /im wavebands. 

3.3. Spitzer Flux Calibration and Color Corrections 

The dominant process that affects the IRAC calibra- 
tion is the location-dependent photometric flat-field cor- 
rections for compact sources, caused by geometric array 
distortions and by the change in the effective filter band- 
pass a s a function of the incidence angles of the radi- 
ation l|Quiiada et all 12004ft . This effect can amount to 
a 10% peak-to-peak calibration error across the array. 
The ultimate absolute point source calibr ation attain- 
able by IRAC is « 2% (iReach et al.ll2005|) . Because of 
the dithering and oversampling used to obtain our im- 
ages, we believe their relative calibration may be closer 
to this latter limit. The complex backgrounds and the 
level of confusion in dense stellar regions of the SMC, 
however, are probably the limiting factors in our photo- 
metric accuracy. 

The MIPS flux densities have an estimated absolute 
uncert ainty of 10% and 20%, for 24 and 70 /im, respec- 
tively ijRieke et al.ll2004|) . These uncertainties are domi- 
nated by the uncertainty in the fluxes of the calibration 
sources. 

To convert our photometric fluxes to the Vega mag- 
nitud system we use the zero points determined by 
M. Cohen (priv. comm.), which are consistent with 
Reach et al. ( 2005). They are listed in Table Q] 

The definition of monochromatic flux density used by 
IRAC, which is common to IRAS and MSX, assumes an 
underlying flux density F v cx i/ -1 . The MIPS flux scale, 
however, assumes a Rayleigh- Jeans underlying flux den- 
sity F v oc v 2 . To convert the MIPS scale into the IRAC 
scale it is necessary to apply factors of 1.041, 1.089, and 
1.043 to the MIPS fluxes measured in the 24, 70, and 
160 /im wavebands respectively. In the magnitude scale, 
these color corrections correspond to -0.043, -0.094, and 
-0.046 mag respectively. Note that throughout this pa- 
per we use the flux definitions native to each instrument. 
Thus these corrections are not applied. 
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4. POINT SOURCE PHOTOMETRY 

Point source photometry was carried out on the mo- 
saic images using the Astronomical Point Source Ex- 
tractor (APEX) software in the MOPEX package (ver- 
sion 050905) provided by the SSC. Using the routine 
apex_lf rame, we detected bright sources in the images 
with a default point response function (PRE) generated 
from in-orbit checkout observations. We selected ~ 30 
of the most isolated of these stars, covering a range of 
brightnesses and positions, and built a new PRF for each 
band with the prf ^response routine. We then fitted the 
new PRF to the data and calculated flux densities for all 
sources detected at a signal-to-noise ratio greater than 
5 (again with apex_l frame). Aperture correction fac- 
tors of 1.12, 1.01, 0.91, 0.92, 0.98, and 1.00 were applied 
to the 3.6, 4.5, 5.8, 8.0, 24, and 70 /im PRF fitting pho- 
tometry respectively, to bring the measured flux densities 
into agreement with very large aperture measurements of 
chosen bright stars. Finally, we created a point source- 
subtracted image with apex_qa to assess the results of 
the PRF fitting and to examine the extended emission 
that remained. 

We can use the source count histograms to estimate 
the fluxes at which incompleteness in the photometry 
becomes significant. Figure |21 shows the histograms of 
the counts for a given source flux in each Spitzer wave- 
band. Assuming that the intrinsic source counts approx- 
imately follow a power law, we use the peak of the his- 
togram to find the flux below which the photometry is 
missing a large number of sources. Because of the effects 
of extended backgrounds and source confusion, the ac- 
tual level of completeness of the photometry depends on 
the location in the image (see i)4.2f> . Nevertheless, Fig- 
ure |3 shows that below 31, 24, 80, 80, 400, and 40000 
/iJy the point source photometry for the 3.6 through 70 
/im wavebands is clearly incomplete. These limits are 
very close to 10 times the 1-a noise in the images, except 
for the 3.6 and 4.5 /jm wavebands where that criterion 
would lead us to expect completeness limits of 7 and 11 
/zJy respectively. The diminished effective "sensitivity" 
in these wavebands is due to the effects of crowding and 
confusion, which are most severe at the shortest wave- 
lengths. 

As in other Spitzer imaging observations, we have 
found that the largest source of false stellar detections 
are t he diffraction artifacts associated with bright stars 
('e.g. JBabbedge et al.ll2006|) . Wherever possible we have 
pruned them from the final catalog, either by hand or 
by a two-step process in the photometry, where we use a 
different PSF for bright stars. In the end, however, the 
catalog suffers from some artifact contamination that is 
difficult to estimate. These artificial sources are usually 
clustered near bright stars, and are generally detected in 
just one waveband. 

4.1. Comparison to Previous Infrared Missions 

In this section we compare the fluxes of different point 
sources in existing catalogs with our Spitzer photometry. 
The SMC was imaged by previous infrared satellites, in- 
cluding IRAS, MSX, and ISO. There are numerous cat- 
aloged sources in the SMC, but the lower angular res- 
olution of most previous missions (especially IRAS) can 
make it difficult to confidently identify them with sources 



seen in our data or at other wavelengths. Also, the dif- 
ferent instrument bandpasses should be considered for 
a fully meaningful comparison of fluxes. In the follow- 
ing discussion we concentrate on the samples for which 
source confusion and bandpass issues are minimized. 

4.1.1. MIPS Comparison 

Usin g the IRAS images, iCroenewegen fc Blommaertl 
(1998) located 30 likely AGB stars in the SMC. Approx- 
imately one third of these sources are carbon stars, one 
third are oxygen-rich AGB stars, and the remainder are 
presumed to be other types of evolved stars. The IRAS 
25 fim bandpass is fairly similar to the MIPS 24 /im band- 
pass, allowing us to compare the flux densities we observe 
with Spitzer to the earlier I RAS observations. Nineteen 
of the stars discussed by ICroenewegen fc Blommaertl 
( 1998) lie within the area covered by our MIPS 24 /im 
map, and 17 have good fluxes in the IRAS and MIPS 
data. 

We find that on average, the IRAS flux densities are 
37% higher than the MIPS flux densities (23% higher if 
we consider the median rather than the mean), and for a 
few sources the difference is as much as a factor of three. 
Only two of the flux densities agree within the quoted un - 
certainties g i ven by iGroenewegen fc Blommaertl l)1998l) . 
iTrams et all |l999) also found a similar discrepancy in 
comparing IRAS and ISO photometry for sources fainter 
than 0.2 Jy. So me of these differen ces are certainly 
due to variability ( Werner et alJ fl98D^ . which in dust- 
shrouded AGB stars can amount up to a magnit ude at 
mid-infrared (MIR) w avelengths (|Le Bertrdll99l Il993t 
Ivan Loon et al.lH998|) . Variability may introduce a sys- 
tematic offset in the sense of one data set consistently 
yielding brighter measurements if that data has consid- 
erably lower sensitivities and the sources in question are 
close to the detection limit (in other words, they are only 
detected by IRAS in their high state). Confusion within 
the large IRAS beam may also play a role, although 
sources that are clearly isolated in the Spitzer observa- 
tions show similar discrepancies between their IRAS and 
MIPS fluxes, suggesting that confusion is not the main 
cause of the disagreement. 

The MSX Point Source Catalog l|Egan et alJ2003[) con- 
tains 9 very bright (> 1 Jy) 21 /im sources in the SMC. 
Eight of these lie within our survey area, and seven are 
unsaturated. Two of these split into multiple sources by 
MOPEX. Comparison of the MSX 21 /im flux densities 
with the MIPS 24 /im flux densities for the remaining five 
sources shows that one of the sources agrees within l-cr, 
and the remaining four are somewhat brighter (13—45 %) 
in the MIPS data. Overall, the average flux ratio be- 
tween MSX and MIPS is 0.85, and the median is 0.88. 
Given the substantial difference between the MIPS and 
MSX bandpasses, and again allowing for variability, it 
seems that the calibration of the 24 /im Spitzer data is 
consistent with that of MSX. 

4.1.2. IRAC Comparison 

At 8 /im there are about 120 MSX sources in our field. 
After matching them to the IRAC catalog and remov- 
ing the sources that are flagged by MSX as being vari- 
able, confused, or low quality, 53 sources remain. These 
sources are on average ~ 10% brighter in our data than 
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Fig. 2. — Source count histograms for the different Spitzer wavebands in black (bin size is 0.2 mag). The bottom axis scale, stellar 
magnitudes, is shared between panels. The top axis scale, stellar flux, is different in each panel due to the different magnitude zero points. 
The photometric completeness thresholds estimated based on 10 times the sensitivities of the images are indicated by the dashed vertical 
lines: 7, 11, 80, 65, 430, and 50000 /ijy for the 3.6 to 70 (im wavebands respectively. The solid lines indicate the peak of the source count 
histograms, below which the photometry is clearly incomplete: 31, 24, 80, 80, 400, and 40000 /ijy respectively. The completeness of the 
photometry in the 3.6 and 4.5 fim wavebands, where the source densities are higher and confusion is more severe, is consequently worse 
th an ex pected based on the sensitivity of the images alone. The white histograms show the source counts for the classes of sources identified 
in H5.5I A large fraction (~ 50%) of the brightest sources in each waveband fall into one of the identified classes. 



in MSX. This effect is entirely attributable to the fact 
that the IRAC 8 /im bandpass is slightly bluer than that 
of MSX {\irac ~ 7.74 /im compared to \ M sx ~ 8.28 
/im). A stellar blackbody in the Ray leigh- Jeans limit 
(valid for photospheric temperatures T p h > 1000 K) will 
be approximately (8.28/7.74) 2 w 1.14 times brighter in 
the IRAC waveband. Thus, the 8 /im IRAC point source 
calibration is consistent with MSX. 

There is also some overlap between the ISO Magel- 
lanic Clouds Mini-Survey (Loup et a!., in preparation) 
and the area mapped by our survey. The ISO data were 
taken in the LW2 band, which is centered at 6.7 /im 
and has a FWHM of 3.5 /im. Unfortunately, this band- 
pass does not match any of the IRAC bands very closely, 
but the average of IRAC 5.8 /im and 8.0 /im flux den- 
sities should provide a reasonable approximation of the 
LW2 flux density. We find 114 sources in the ISO cat- 
alog that have both 5.8 /im and 8.0 /im fluxes, have no 
close neighbors, are not confused, and have an ISO S/N 
of at least 10. After discarding the 10 reddest sources 



([5.8] — [8.0] > 0.51), and using the mean of the 5.8 /im 
and 8.0 /im bands, we find that the IRAC flux densities 
are ~5% larger than the ISO measurements with a scat- 
ter of ~33%. Given the uncertainties involved in com- 
paring the different bandpasses, the calibration of our 
IRAC data appears compatible with that of ISO. 

4.2. Photometric Completeness 

We used artificial star tests to estimate the complete- 
ness of the photometry We added scaled versions of the 
measured PRF at random positions in each image and 
re-ran the MOPEX photometry to see how many of the 
artificial sources were recovered. Because of the large 
variations in both the surface density of point sources 
and the brightness of extended emission regions across 
the SMC, the completeness is a complex function of both 
flux and position. In order to determine the approxi- 
mate size of these effects, we carried out completeness 
tests in four separate fields that span the range of sur- 
face densities found in the data. Field 1 is centered at 
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(0:2000, £2000) = (00 h 47 m 53!5, -72°22'46") in the south- 
west part of the Bar, with extremely high stellar den- 
sity and significant extended emission. Field 2 is at 
(o 20 oo,<$2ooo) = (00 h 54 m 33 i ;7, -72°34'23") near the mid- 
dle of the Bar with a high stellar density. Field 3 is at 
(o 20 oo,<W)) = (01 h 01 m 54?l,-72°00'27") near the top 
of the Bar with a moderate stellar density and some ex- 
tended emission. Finally, Field 4 is well out in the Wing 
at (a 20 oo,<520oo) = (01% m 49?2, -73°25'08"), in a low 
density area. Each field is 768 pixels on a side, covering 
an area of 104.9 arcmin 2 in the IRAC bands and 1065.4 
arcmin 2 in the MIPS 24 /im mosaic. 

We summarize in Table El the results of the artificial 
star tests. The worst crowding occurs at 3.6 /im, and in 
the highest density regions the photometry is only 90% 
complete to [3.6] = 16 (110.3 /iJy). in the lowest den- 
sity regions the 90% completeness limit is fainter than 
[3.6] = 17. At longer wavelengths, the crowding declines 
dramatically as the stellar photospheres that comprise a 
large majority of the detected sources fall off as A -2 and 
the detector sensitivity rapidly decreases. In the MIPS 
data, incompleteness is largely caused by the presence of 
bright dust emission rather than point source crowding 
(as there are no very crowded regions). Thus, the com- 
pleteness results in the four fields scale as the percentage 
of the field that is covered by bright extended emission. 

4.3. Catalog Compilation 

We searched for sources in each waveband using APEX 
and the procedure outlined in 2] The catalogs for the 
individual wavebands were then merged using positional 
tolerances of 1" for IRAC, and 3" and 6" for MIPS at 
24 and 70 /im respectively. Catalog assembly proceeded 
from the shortest to the longest wavelength; the closest 



source within the aforementioned positional tolerances 
was associated with a pre-existing source, or a new en- 
try was created if no cataloged source existed within the 
tolerances. We derived the source position in the final 
merged catalog of all the Spitzer bands using the signal- 
to-noise weighted average of the detections in the IRAC 
bands, or adopting the the 24 or 70 /im positions (in that 
order) if no IRAC detections were available. 

We proceeded to cro ss-identify the Spitz er sources with 
the OGLE II catalog <|Udalski et al.lll99^ and the SMC 
catalog of the Magellanic C loud Photometric Survey 
(MCPS, iZaritskv et all 12002ft in the o ptical, and with 
the 2 MASS catalog in the near infrared iSkrutski e et alJ 
I2006|h The fluxes report ed by OGLE II were preferred to 
those of the MCPS fsee IZaritskv et al.ll2002() . The posi- 
tional tolerance used for cross-identification was 1". The 
final catalog compiled has over 400,000 point sources, 
and is available electronically at the project URL 1 . To 
produce a uniform catalog we expressed all fluxes in /iJy, 
using the following photometric zero points for BV IJHK: 
4060 3 723, 2459, 1594, 1024, 666.8 Jy JCohen et alJ 
l2003albf) . For the Spitzer wavebands we use the zero 
points listed in Tabled 

5. DISCUSSION 

The S 3 MC data provide a rich resource for studying 
many aspects of the interstellar medium (ISM) and star 
formation activity in the SMC. The 3.6 and 4.5 /im wave- 
bands are mostly sensitive to stellar photospheres and 
very hot circumstellar dust, with some contribution from 
bound- free transitions, free-free radiation, and Brackett 
a recombination emission in H II regions. Although the 
4.5 /im waveband is free from AIB contribution, the 3.29 
fim emission feature (thought to be due to C-H bond 
stretching in PAHs) contributes to the 3.6 /im intensity 
near molecular clouds. In cool stars the 4.5 /im waveband 
includes the CO fundamental bandhead, which will ap- 
pear in absorption. Both the 5.8 and 8.0 /zm wavebands 
show molecular material in the ISM and circumstellar en- 
velopes, as well as increasingly faint stellar photospheres. 
They encompass the 6.2 /im emission feature and the 
very bright emission complex at 7 — 9 /im thought to be 
dominated by C-C stretching modes (with some contri- 
bution from in-plane C-H bending) of the bonds in PAHs 
and very small carbonaceous dust grains. The 24 /im 
waveband is sensitive to continuum emission from PAHs 
and VSGs, which are stochastically heated to tempera- 
tures of ~ 150 K by the interstellar radiation field (ISRF) 
or by the central star in warm circumstellar envelopes. 
Thus, the extended 24 /im emission is bright in dusty re- 
gions near massive stars and stellar clusters. In material 
that is shocked or exposed to very hard radiation fields 
associated with peculiar objects, the 24 /im waveband 
may also have a contribution from spectral lines, such as 
[O IV] at 25.89 /im. Emission in the 160 /im waveband 
arises mostly from large dust grains in the ISM in equilib- 
rium with the radiation field, which typically heats them 
to temperatures of several tens of degrees Kelvin. The 70 
/im waveband, finally, is a combination of VSG and large 
grain emission, in proportions that depend on the rela- 
tive a bundances of the grain populations (jDesert et alJ 
1990). 

1 http : //celestial . berkeley . edu/ spitzer 
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5.1. The SMC in the Mid and Far Infrared 

The Spitzer images of the SMC reveal an astounding 
degree of detail. Figure[3]is a color composite of the 8.0, 
4.5, and 3.6 /ini wavebands corresponding to the red, 
green, and blue colors. The regions displaying bright red 
emission are associated with star formation and known 
molecular clouds: N 66 (NGC 346), the region around 
N 76, N 83/N 84, and the south-west tip of the Bar are 
all prominent at 8.0 /im. The south-west region is also 
the region of highest stellar density, as evidenced by the 
bright diffuse blue background made up of sources be- 
yond Spitzer 's confusion limit. Note that because the 
3.6/5.8 /xm and the 4.5/8.0 /im wavebands of IRAC oc- 
cupy different places in the Spitzer focal plane, the areas 
mapped by each pair of channels do not completely over- 
lap. 

Close inspection of Figure reveals that a dramatic 
change of the colors of the MIR extended emission takes 
place within star forming regions. Indeed, Figurc^Jshows 
clearly that the bright extended 8.0 /xm emission fre- 
quently surrounds a well defined cavity around one or 
more stars, inside which the emission is bluer. This phe- 
nomenon is likely caused by the destruction of the carri- 
ers of the 8.0 /mi emission (nominally PAHs), which are 
obliterated by the intense radiation field present in H II 
regions. The transition between the regions of intense 
and faint 8 /xm emission is often sharp, probably corre- 
sponding to the edges of dense molecular material where 
shielding protects PAHs from photodestruction by UV 
radiation. 

What percentage of the flux in the different wavebands 
is due to point sources? Moreover, to help interpret the 
unresolved infrared observations from distant metal-poor 
galaxies it is interesting to estimate the contribution to 
the global infrared flux arising from various types of point 
sources. TableOllists the fluxes measured over the entire 
S 3 MC images and compares it to the fluxes in the differ- 
ent classes of point sources discussed in *15.5I Considering 
the level of completeness in the phometric catalog in the 
areas of high confusion (e.g., the SW Bar), it is clear 
that most of the flux in the 3.6, 4.5, and 5.8 /xm wave- 
bands is due to point sources. At 8 /xm, however, there 
is a noticeable contribution from extended emission that 
becomes rapidly dominant at the longer wavelengths of 
MIPS. Among the different types of sources discussed in 
M5.5I the carbon stars stand out prominently in the IRAC 
wavebands, while the YSOs are the largest fraction of 
resolved flux for MIPS. The red supergiants and oxygen- 
rich AGB stars, and the dusty early-B stars, although 
clearly detected in the analysis, are not prominent con- 
tributors to any of the Spitzer bands. 

5.2. PAHs in the SMC 

Among the recent Spitzer results is the disappear- 
anc e of PAH emission in galaxies of low metallic- 
ity llEngelbracht et al.ll2005t iHogg et al.ll2005l IWu et all 
120061 IQ'Halloran et a.U 12006^ The F 8 /F 24 flux ra- 
tio in galaxies changes dramatically at a metallicity of 
12 + log(0/H) ~ 8 (the metallicity of the SMC), with low 
metallicity objects exhibiting very low ratios (Fg/F24~ 
0.08) compared to those seen in solar-metallicity ob- 
jects (Fg/F24^ 0.7). This decrease in the 8 /xm emis- 
sion relative to 24 /xm is frequently interpreted as the 



enhanced destruction or reduced formation of PAHs 
in th e se environments llPlante fc Sauvagell2002t iMaddenl 
12001 iMadden et alJl2006|) . PAH destruction by itself, 
however, is probably not enough to account for the ob- 
served range in the F8/F24 ratio, which may require also 
an en hanced hot dust compon ent in low-metallicity sys- 
tems l)Engelbracht et aT1l2005() . Figures and 0] clearly 
show that there are large variations in the brightness of 
the 8 /xm emission depending on the local conditions. 
These images suggest on the sca le of the entire SMC an 
important conclusion reached bv lReach et al.l 1)2000]) us- 
ing observations of a small, quiescent region: PAHs exist 
even in the low-metallicity interstellar medium when the 
local environment allows it. 

We further explore the variation in the relative inten- 
sity of the 8 /im waveband by producing a ratio map of 
the SW region of the Bar (Figure [SJ. In order to pro- 
duce this image we convolved the 8 /xm image to a 6" 
resolution, subtracted the foregrounds from each image 
(estimated from the modes of the respective histograms 
to be 5.09 and 21.90 MJy sr" 1 for the 8 and 24 /tm 
images respectively), and median filtered the resulting 
images to remove most point sources. It is apparent that 
the Fg/F24 ratio takes a large range of values, with an 
overall gradient increasing toward the southwest. This 
gradient has the opposite sense to the gradient in the 
brightness of the 24 /im emission, an indicator of mas- 
sive star formation and intense ISRF. Extended regions 
of the SW Bar exhibit F 8 /F 24 > 0.7 typical of solar- 
metallicity objects, while other areas generally associ- 
ated with bright 24 /xm continuum and massive star for- 
mation have F§/F 2 4 < 0.1. Because the re are no larg e 
variations of metallicity across the SMC (Dufour 1984), 
it appears that the PAH abundance (at least as traced by 
the Fg /F24 ratio) is regulated primarily by the local ISRF 
rather than metallicity. This does not imply that metal- 
licity plays no role on setting the abundance of PAHs on 
galactic scales. In the context of global measurements of 
galaxies, where the ISRF is a function of metallicity, this 
observation suggests that photodestruction of PAHs, not 
reduced formation, may be primarily responsible for the 
range of values measured for the F§/F24 ratio. 

Figure shows a quantification of the relationship be- 
tween this ratio and the dust FIR color F24/F70. To 
obtain this plot we convolved the 8 and 24 /xm images 
of the SW Bar region to the angular resolution of the 
70 /*m data, resampled the median-filtered maps to ob- 
tain approximately one pixel per resolution element, and 
removed the corresponding foregrounds from the indi- 
vidual images (in addition to the 8 and 24 /xm fore- 
grounds mentioned above, we estimate the 70 /xm fore- 
ground to be ~ 6.5 MJy sr" 1 ). We considered only 
data with F 8 > 0.2, F 24 > 0.2, and F 70 > 6 MJy sr" 1 
after foreground removal. We computed the density 
of points in the F§/F24 ratio versus the F24/F70 ratio 
plane, and calculated its mean in F24/F70 bins. We 
find that, albeit with large scatter, the trend for the 
Fg/F 2 4 ratio as a function of F24/F70 is described by 

log[F 8 /F 24 ] 0.71og[F 2 4/F 7 o]-1.4. The uncertainties 

in these parameters are considerable and due mostly to 
uncertainties in the foreground substraction. This shows 
that regions with blue dust colors (thus presumably high 
radiation fields) tend to be underluminous in 8 /xm emis- 
sion with respect to 24 /xm, again suggesting a link be- 
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Fig. 3. — Color composite of the S 3 MC mosaic for the IRAC 8.0, 4.5, and 3.6 /im wavebands, represented by the colors red, green, 
and blue. The prominent regions are N 83/N 84 at the SE tip of the Wing, N 76 in the north, N 66/NGC 346 near R.A.= l h 00 m and 
Dec.= —72° 10', and the entire southwest region of the Bar, which hosts several prominent H II regions such as N 22. Because of the 
arrangement of the different IRAC wavebands in the focal plane of Spitzer, the 3.6 fim data is displaced westward from the 4.5 and 8.0 fim 
data. 



TABLE 3 
Resolved Flux Statistics 



Waveband 


Total Flux 1 

(Jy) 


Point sources 

(Jy) 
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(%) 


Carbon stars 
(%) 


RSG/OAGB 
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Dusty B stars 
(%) 
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0.0 


4.5 fira 


115 


71.4 


0.4 


14.4 


3.2 
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17.2 


3.2 
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8.0 /im 
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5.0 


18.6 


2.3 


0.6 


24 fim 
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73.3 


26.6 


2.0 


0.3 


4.4 


70 /im 


9451 


347 


19.7 


0.1 


0.0 


3.4 2 



1 Approximate total flux in the corresponding Spitzer image after foreground subtraction, 
that the mapped area is different in the various wavebands. 
2 From only 6 objects detected. Probably misidentified 70 /im sources. 



Note 



tween the intensity of the local ISRF and the presence 
of PAH emission. The determination of the precise roles 
of PAH destruction and enhanced 24 /im VSG emission 
in setting the F s /F 2 4 ratio will have to await upcoming 



spectroscopic observations. About 55% of the 8 fxm flux 
and 22% of the 24 /im flux in the SW Bar image are as- 
sociated with regions where Fg/F24 > 0.5. The overall 
ratio of the total 8 and 24 /im fluxes in this region of the 
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Fig. 4. — Details of Figure|2]on four star forming regions. The sharp change in the MIR colors in the region surrounding ionizing sources 
is particularly dramatic for N76 and N83/N84. These regions are mostly devoid of 8 £im emission. 



SMC is F 8 /F 24 « 0.36. We defer the detailed quantita- 
tive study of the extended emission of the SMC in the 
FIR and MIR to a series of subsequent papers (Leroy et 
al., in prep.; Sandstrom et al., in prep.; Bolatto et al., in 
prep.). 

5.3. Infrared Point Sources and Their Optical 
Counterparts 

Perhaps one of the most interesting matters to pursue 
using the S 3 MC point source catalog is to investigate 
the classes of sources detected by Spitzer in its different 
wavebands. The cross identification against optical data 
is very valuable in this regard, as the techniques of source 
classification based on the MIR colors and magnitudes 



alone are currently in their infancy. The S 3 MC catalog 
provides unique deep photometric data on a large sample 
of stars located at the same distance, with similar metal- 
licities, and low levels of foreground and background con- 
fusion. We expect that one of the contributions of this 
dataset will be to provide a basis for testing different 
source-identification schemes. 

Figure [7] shows the optical color-magnitude diagrams 
(CMDs) for sources detected in each Spitzer waveband 
at a signal-to-noise greater than 10. Several popula- 
tions are immediately apparent: 1) the main sequence 
for stellar types earlier than AO, at approximately zero 
V — I color; 2) the blue-loop helium core-burning stars, 
running parallel to the main sequence; 3) the red giant 
branch and the red clump, where most of the Spitzer 
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Fig. 5. — The 24 (im image of the SW region of the SMC Bar at 6" resolution, and the corresponding Fg/F24 ratio, an estimator of 
PAH abundance. The units of the 24 (im intensity scale are MJy sr _1 , and the display is saturated in the bright regions to emphasize the 
low level emission. Regions with surface brightness below 0.2 MJy sr _1 in either band are mask ed out of the ratio. There is an overall 
anticorrelation between bright 24 /an emission and large Fs/F24 ratios. The region studied by [Reach ct al. (200(f), where PAHs were 
detected by ISO toward the molecular cloud SMCB1#1, is outlined by a rectangle (white in the 24 fim panel, black in the ratio panel). 
The area of high ratio within that region looks essentially identical to the 11.3 fj,m AIB map produced with ISO's CVF, suggesting that 
the Fg /F24 ratio is an excellent diagnostic of PAH abundance. 




Fig. 6.— The F 8 /F 2 4 ratio in the SW region of the SMC Bar 
as a function of the F24/F70 ratio. The black squares indicate 
the averages of the F8/P24 ratio in bins of F24/F70. There is 
a clear anticorrelation between the Fg/F24 and F 2 4/F7o ratios, 
showing that regions with very blue dust FIR colors are most likely 
underabundant in PAHs. The dashed line shows a power-law fit to 
the data, log[F 8 /F 24 ] ~ -0.71og[F 2 4/F 7 o] - 1.4 



sources are found, at V — I ~ 1; 4) the asymptotic gi- 
ant branch (AGB) and the red supergiants (RSG); 5) 
the approximately horizontal branch corresponding to 
evolved carbon stars, at My ~ —2 and V — I > 2; 
and, 6) a fuzz of sources at low luminosities that consist 
of misidentifications, stars in the Milky Way halo and 
thick disk, unresolved background galaxies, and highly 



reddened sources. The location of the blue-loop stars 
is both theoretically and observationally linked to the 
metallicity of the intermediate-mass stars that populate 
this feature; i.e., theoretical isochrones show a metal- 
licity dep endence on the location of the blue loop i n 
the CMD lIGirardi et alJ2000tlLei eune fc Schaereil200l . 
and the observations of different metallicity dwarf irreg- 
ular galaxies confirm that the location of the blue-loop 
stars varies with metallicity as predi cted by isochrones 
Ce.g.. lDohm-Palmer et al1fl997l [1998^ . In the SMC, the 
blue loop stars are well separated from the main-sequence 
as expected for their metallicities. The foreground plume 
is nearly vertical at V — I = 0.7 whereas the RSG stars 
occupy a locus that is slightly slanted atV — J ~ 1.1 — 1.2 
and My ~ — 3 to almost —6. The foreground plume 
has been seen in other studies of th e stellar populations 
of Lo cal Group dwarf galaxies (e.g. IMcConnachie et alJ 
2005). 

As the sensitivity of our observations, indicated by the 
gray bar on the main sequence, diminishes for the longer 
wavebands, we progressively lose the fainter sources. The 
sensitivity across the CMD, however, is clearly dependent 
on the V — I color and follows a trend, illustrated by the 
dashed gray lines in Figure [7] that is approximately par- 
allel to the reddening vector illustrated in the upper right 
corner of each panel. This reddening vector corresponds 
to Ay/E(V — I)~ 1.85, derived using SMC measure- 
ments and assuming a Galactic E(V — I)/F,(V — J)tv 0.72 
l|Bouchet et all 1 9851 iRieke fc Lebofskvll985|) . The color 
dependence of the sensitivity suggests that many sources 
in the highly reddened region of the CMD could be in- 
trinsically brighter stars located behind a few magnitudes 
of extinction. The spatial distribution of these reddened 
stars, however, is not strongly clustered. Such clustering 
would be expected if they were located behind localized 
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Fig. 7. — Optical color-magnitude diagrams of the stars detected in the different bands of Spitzer for which optical conterparts were 
found in the OGLE II and the MCPS catalogs. The overlaid gray contours indicate density of sources. The assumed distance modulus 
is DM = 18.93, and the photometry has been corrected by the foreground reddening and extinction toward the SMC (Ay ~ 0.12, 
E(V — I) ~ 0.09). The right side scale in the last panel correspond to apparent magnitudes. The gray bars indicate the expected 10-cr 
detection thresholds for naked sources on the main sequence in the different Spitzer wavebands, given the IRAC point source sensitivities 
(c.f., Table HI and the known color corrections. The vector in t he upper right corner of the first panel indicates the extinction correction 
for one magnitude of reddening, assuming Ay /E(V — I)= 1.85 IBouchet et al.lll985l) . The dashed gray line in the 5.8 and 8.0 /an panels 
is parallel to the extinction vector and indicates the boundary over which the density of detected sources drastically increases. No naked 
photospheres should be detectable in the MIPS data: detection at 24 and 70 fim requires extensive reprocessing of the photospheric 
radiation by dust. All the sources in the 70 fim diagram are consistent with misidentifications (their distribution in the CMD plane is 
the same as for the 3.6 or 4.5 fim sources). The rightmost panel sketches the loeii of the different branches with their identification: MS 
(main sequence), BL (blue-loop, helium core-burning stars), FG (foreground), RSG (red super giants), OAGB (oxygen-rich asymptotic 
giant branch), C* (carbon stars), RGB (red giant branch). 



regions of extinction, such as molecular clouds. Clearly, 
if the reddening of these stars is caused by dust, this dust 
must be closely associated with the sources themselves. 
Most of these sources, however, do not exhibit a strong 
infrared excess, which makes this explanation unlikely. 
Alternatively, these sources may be K and M dwarfs in 
the disk and halo of our own Galaxy, at distances of 
less than 2 kpc. Their infrared colors (J — H ~ 0.6, 
H — K ~ 0.2) support this hypothesis. The slope in the 
sensitivity would then be related to the slope of the lower 
main sequence, weighted by the distance distribution of 
the sources. In addition to these foreground interlopers, 
there is contamination by background, unresolved galax- 
ies, which can be seen throughout the SMC at infrared 
wavelengths. We discuss the importance of these con- 
taminants in i)5.5.1l 

We do not expect to detect stellar photospheres in the 
MIPS wavebands. Indeed, the drops in detector sensi- 
tivity and in stellar flux (which approximately follows 
a A~ 2 Rayleigh- Jeans law), combine to lower our detec- 
tion limit by ~ 4.4 mag at 24 /xm with respect to the 
8.0 /im IRAC waveband. Thus, only the brightest red 
supergiants and AGB stars are directly detectable at 24 
/xm. The other sources found in the MIPS images must 
have a substantial infrared excess from dusty envelopes or 
emission lines, or correspond to misidentifications. The 
latter is probably the case with the optical counterparts 
of most sources at 70 /*m, which are distributed in the 
CMD plane very much like a randomly drawn subsample 
of the population detected at 3.6 or 4.5 /xm. 

Besides the bright end of the AGB and the RSG stars, 



there are three other distinct populations of sources 
clearly detected at 24 /im. The bright end of the 
main sequence shows an overdensity of sources with 
V — I ~ —0.15 and absolute V magnitudes in the range 
—2.5 > V > —4.5, corresponding to stellar types B3 
to 09. The highly reddened side of the diagram shows 
the branch corresponding to evolved carbon stars, with 
an overdensity of sources connecting it to the red super- 
giants. Finally, the spread of sources at the bottom of 
the diagram, which persists also throughout the IRAC 
bands, is not entirely consistent with a randomly drawn 
subsample of the detections at 3.6 or 4.5 /xm. Although 
there is probably strong contamination, if the bottom of 
the CMD were entirely composed of misidentified sources 
we would expect it to look much more like the 70 /xm dia- 
gram, with a source density following closely the contours 
of the CMDs for the lowest IRAC wavebands and featur- 
ing a prominent red clump. It does not, thus a sizable 
fraction of these sources are indeed bright MIR sources 
with very faint optical counterparts. It turns out that 
many of them have very red MIR colors, and are consis- 
tent with background galaxies and foreground M dwarfs. 
We will investigate this matter further in H5.5I 

The fact that the carbon star branch is prominent at 
24 /xm suggests that it may be possible to use Spitzer ob- 
servations to identify these sources, which are common 
contaminants in extragalactic chemical evolution studies. 
Carbon stars are difficult to analyze due to their plethora 
of carbon molecular lines, which affect the continuum 
placement in the analysis of the spectra. This problem is 
even more accute in low-metallicity galaxies, where car- 



14 



bon stars are more abundant llBlanco e t al.lll978t 119801) 
and have stronge r C2H2 bands ijvan Loon et alJ 12003: 
iSloan et"al] 12006^ . Recent work on dwarf spheroidal 
galaxies, for example, has been b ased on detailed spec- 
tral analyses of RG B stars (e.g., iShetrone et alJ 120031: 
iTolstov et all 1200.1 iSadakane et al J 12004) . The selec- 
tion of the target stars, however, is not always obvious 
and suffers from the inclusion of carbon stars. In fact, 
because the studied galaxies are far away the stars se- 
lected tend to be at the tip of the RGB, where there is 
the greatest overlap and confusion potential with carbon 
stars. In the aforementioned studies these contaminants 
amounted to 25% of the original sample, and weeding 
them out required a large time investment in an 8m class 
telescope. With infrared colors, tar get selection can b e 
significantly improved. As shown bv lSloan et alJ {2006), 
and also seen here, the carbon stars easily separate out 
from RGB stars at infrared wavelengths, e.g., 3.6, 4.5, 
5.8, 8.0, and even 24 /mi. Analysis of the old RGB stars 
in the dwarf irregular galaxies may show similarities to 
the metal-poor halo of our Galaxy, and help to explain 
the formation of our Galaxy, unlike the chemistry of the 
metal-poor stars in the nearby dwarf spheroidal galaxies 
which do not show chemic al similarities with metal-poor 
stars in our Galaxy at all ijVenn et alJl2004j) . 

5.4. Mid and Far IR Color- Magnitude Diagrams 

Figure |H1 shows the color-magnitude diagrams in dif- 
ferent combinations of the Spitzer wavebands. The pro- 
gression from left to right shows the absolute magnitude 
in one band plotted against the color with respect to the 
next Spitzer waveband. The rightmost panel shows the 
8 /im magnitude versus th e [3.6] — [8.0] color, a CMD 
modeled bv lWhitnev et alJ (|2003|) . 

Several features are apparent in these plots. Since 
most stars emit like blackbodies in the infrared, they 
pile on a vertical line near zero color constituting a MIR 
"main sequence" . At longer wavelengths there is an in- 
creasing number of objects that appear as a red plume 
of sources, which becomes progressively more detached 
from this main sequence. In several of these CMDs there 
is also a number of sources in a branch off the main se- 
quence toward redder colors at high luminosities. At the 
longer wavebands the progressive loss of sensitivity re- 
sults in a diminishing number of sources detected on the 
main sequence, until in the 24 /mi vs. [24] — [70] CMD 
most sources have a color redder than 3, implying a ra- 
tio of fluxes F70/F24 > 1.8 and a spectral energy density 
rapidly rising toward the far infrared. 

5.5. Stellar Populations 

To gain insight into the nature of the different stellar 
populations detected, we used color and flux selection 
criteria in the CMDs to pull out particular subsets of 
sources (Figures El and ■ We defined four groups of 
sources, which we discuss below. Sources in groups 1 to 
3 are a substantial fraction of the brightest sources seen 
in our Spitzer data (> 50%, FigureEJ). 

1. Sources with [5.8] - [8.0] > 1.2, and -10 < M 5 . 8 < 
—6. We also require that these sources be detected 
with S/N > 10 at these wavebands and one neigh- 
boring waveband (either 4.5 or 24 /mi). About 280 



sources fall into this category. They are represen- 
tative of the highly reddened plume present in the 
5.8 /mi CMD. 

2. Sources with 0.16 < [5.8] - [8.0] < 0.6 and -12 < 
M 5 . s < -7.5 and S/N > 10. About 660 sources 
fall into this class. These correspond to the branch 
that splits off the upper main sequence in the MIR. 
These sources are identified by their optical CMD 
as carbon stars. 

3. Sources detected with S/N > 10 at 24 /mi and 
with optical colors 1.2 < V — I < 1.9 and absolute 
magnitude —6 < My < —5. About 130 sources 
belong to this category. These represent the red 
supergiants and the bright end of the OAGB (i.e., 
oxygen-rich evolved stars) detected at 24 /mi. 

4. Sources detected with S/N > 10 at 24 /mi and at 
3.6 or 4.5 /mi, with optical colors —0.4 < V — I < 
and absolute magnitude —5 < My < —2.5. These 
sources correspond to the population of early-B 
stars (B3 to 09) detected at 24 /mi; about 190 
(5%) of the stars that fulfill these optical color 
criteria have 24 /im counterparts. Inspection of 
the individual SEDs to cull sources with poor or 
likely misidentified photometry reduces the sample 
to ~ 120 sources. 

These populations have been highlighted in the optical 
CMDs in Figure El This plot shows that the stars with 
large 8.0 /im excess selected in the MIR CMDs (group 1 
defined above) have only very faint optical counterparts, 
as they are almost all located in the spread of sources at 
the bottom of the diagram. Many of these counterparts 
may in fact be misidentifications. The bright stars with 
small 8.0 /mi excess (group 2) are clearly identified with 
the branch of carbon stars in the optical. Bright car- 
bo n stars have been studied using Spitzer spectroscopy 
bv ISloan et alJ (|2006f) . They find that oxygen-rich and 
carbon-rich AGB can be split into two sequences us- 
ing J — K and [8] — [24] colors, with the oxygen-rich 
and carbon-rich sources forming two sequences around 
1< J - K < 2 and [8] - [24] - 1 - 3, and J - K > 2 an d 
[8] - [24] - 0.8 respectively (see also lBlum et~aTll2006|) . 
We clearly see both sequences in Figure ll2Ti. The car- 
bon stars, as well as the red supergiant stars selected in 
group 3, are not preferentially distributed along the SMC 
Bar but show a much smoother spatial distribution (Fig- 
ure lllbx). In particular, the distribution of the MIR- 
selected carbon stars is compatible with the spheroidal 
shape found by other studies of the older SMC stella r 
populations fe.g.. lCioni et al.l2000l:lZaritskv et alJ 2000). 

Identification of these populations in the CMDs in Fig- 
ure [TO| shows some of their MIR characteristics. The red 
supergiants, for example, are also most of the bright- 
est stars in the MIR. Many of the rest of the brightest 
stars are carbon stars, which form a distinct branch off 
the main sequence not only in the [5.8] — [8.0] color, but 
also in the K - [3.6] CMD. The B stars with 24 /mi 
excess cluster near the bottom of the main sequence in 
most of these CMDs due to their large bolometric correc- 
tions at IR wavelengths, although they tend to be under- 
represented at the longest wavelengths because very few 
are detected at 5.8 or 70 /im. 
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Fig. 8. — Spitzer infrared color-magnitude diagrams of the SMC. The waveband corresponding to the ordinate axis is indicated in each 
panel. The overlaid gray contours indicate densities of sources, and the assumed distance modulus is DM = 18.93. Only sources detected 
with signal-to-noise greater than 10 in the wavebands used in each CMD arc ploted. The right side scale on the last panel indicates apparent 
magnitude. 
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FlG. 9. — Same as FigureUI now with different populations selected by the magnitude and color cuts defined in the main text highlighted. 
The green, red, blue, and magenta symbols correspond to populations number 1 (YSOs), 2 (carbon stars), 3 (RSG and OAGB stars), and 
4 (B stars with 24 fim emission) respectively. Most of the stars in population 1, sources with strong 8.0 fim excess, have only very faint 
optical counterparts. Sources in population 2, which splits off the bright end of the MIR main sequence, are clearly identified with the 
branch of carbon stars. 



5.5.1. A Sample of Young Stellar Objects in the SMC 

The highly reddened 8.0 /jm sources in group 1 consti- 
tute many of the sources detected at both 24 and 70 /im. 
Their spatial distibution fFigureHTk') is highly clustered 
around known molecular clouds and star forming regions, 
suggesting that many of these sources are still embedded 
young stellar objects (YSOs), with a secondary contribu- 
tion from unresolved extragalactic background sources. 
Further support for this hypothesis is lent by the loca- 
tion of these sources in the [8.0] vs. [3.6] — [8.0] CMD, as 
well as in the color-color plots in Figure These com- 
binations of colors have been modeled by W hitney et al.l 



( 2003) for different objects, including YSOs. The highly 
reddened 8.0 fim sources occupy the space in those plots 
that is predicted for late Class to Class II YSOs with 
different geometrical parameters. About 280 candidate 
YSOs fulfill our color-magnitude criteria, and they are 
compiled in Table 4 (published in its entirety in the elec- 
tronic edition of the journal). The color-magnitude se- 
lection criteria used to identify these sources have been 
chosen to minimize the contamination of the sample, and 
limit us to the brightest YSOs in the SMC; clearly there 
is a large population of fainter YSOs detected in this sur- 
vey as well, but it is highly contaminated by background 
galaxies (Figure ITTk). By modeling the SEDs of individ- 




ual sources rather than relying only on color selection it 
is possible to identify a larger samples of YSOs in these 
data. Simon et al. (in prep.) successfully use this tech- 
nique to identify over 200 YSOs with stellar masses as 
low as 1.5 Mq in the N 66 region alone. 

Missing from Figure El are the early Class YSOs, 
which are predicted to have extreme [24] — [70] colors 
> 12. Our sensitivity at 24 /xm would not, for the most 
part, allow us to detect such sources. Thus, some of 
our 70 /jm detections without 24 /im counterparts may 
be early Class YSOs (Figure Ef). There are approx- 
imately 370 highly significant (S/N > 25) 70 /jm de- 
tections without 24 /xm counterparts in our photometry, 
with a total flux F70 ~ 127 Jy (i.e., almost twice the total 
flux in our sample of YSOs in Table l5~5~Tl YSOs as noted 
in Table [3Jl . Because of the limited angular resolution 
at 70 /xm, the image artifacts, and the bright extended 
emission result in catalog reliability and confusion issues, 
however, we defer such a study to a separate paper (Si- 
mon et al., in prep.). 

To quantify the degree of contamination by back- 
ground galaxies and foreground thick disk interlopers we 
have used the SWI RE catalog of the ELAIS-N1 field by 
iSurace et alJ ((2004) . This field is conveniently located at 
the same Galactic latitude as the SMC but in the north- 
ern hemisphere, and at a Galactic longitude I ~ +90° 
comparable to that of the SMC (I ~ —50°). Applying 
exactly the same color and magnitude criteria we used 
to find YSOs in the S 3 MC catalog to the SWIRE catalog 
we find 7 sources in the entire 8.7 square degree ELAIS- 
Nl field l)Lonsdale et al J 120031) . Scaling this number to 
the 5.8 and 8.0 /xm overlap area in the S 3 MC mosaic 
(2.45 square degrees), we expect ~ 20 contaminants in 
our sample of ~ 280 candidates (~ 7%). The magnitude 
and color cuts employed are designed to minimize the de- 
gree of contamination: pushing below M5.8 = —6 rapidly 
increases the number of background galaxies included. 
Reducing the faint magnitude cut to M5.8 = —5 results 
in ~ 110 interlopers out of 533 sources (~ 21%). Moving 



the limit in [5.8] — [8.0] to bluer colors, although more 
benign than including fainter sources, results also in a 
higher rate of contaminants. Using [5.8] — [8.0] = 0.8 (1.0) 
yields 40 (27) interlopers in a total of 386 (340) sources, 
or ~ 10% (8%) contamination. FigurelTTfe shows the spa- 
tial distribution of sources that have the same [5.8] — [8.0] 
colors as our YSO candidate sample but that are fainter 
than M5.8 = —5. Although there is an concentration of 
such sources along the Bar and in known star-forming 
regions, they are much more uniformly distributed than 
our YSO sample, pointing to heavy contamination by 
sources outside the SMC. 

5.5.2. A Population of Dusty Early B Stars in the SMC 

Figure [5] shows that the early B stars detected at 24 
/xm, defined as group 4 above, have an extremely inter- 
esting spectral energy distribution. They gradually dis- 
appear in the 5.8 and 8.0 /xm wavebands, then strongly 
reappear at 24 /xm, and there are virtually no detections 
at 70 /xm. These 24 /xm-bright stars are not conspicuous 
at any other wavelength. Their typical absolute magni- 
tudes in any of the IRAC bands are Mxrac ~ — 4, while 
their mean 24 /xm absolute magnitude is M24 ~ —10. 
Their mean SED (Figure I13fl shows that the typical 
24 /xm excess over the photospheric flux is very large, 
F-z&l Fphoto ~ 200, otherwise we would not have de- 
tected them. Their K — [24] color, an excess indicator 
widely used because of its independence from the stel- 
lar temperature, is extremely large, with a typical value 
K - [24] - 6.5. 

The sharp rise in the SED of these stars at 24 /xm can 
be explained by emission from ~ 0.05 — 1.4 M ffi of dust 
at a temperature of Tdust ~ 150 — 200 K, or perhaps 
by an emission line present in the 24 /xm window. It is 
difficult to imagine what such an emission line could be, 
since the usual candidates such as [O IV] at 25.89 /xm or 
[Fe II] at 24.52 /xm are not prominent in normal stars. If 
dust is the cause of their detection at 24 /xm, however, 
a particular geometry is required because these stars are 
not appreciably reddened in the optical. Their intrinsic 
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Fig. 11. — Distribution throughout the SMC of the different populations of sources selected using the color-magnitude cuts described 
in the main text: (a) YSOs. (b) Carbon stars, (c) Red Supergiants and bright end of the AGB. (d) Early B stars with 24 /an emission, 
(c) Faint and red MIR sources, corresponding to YSOs heavily contaminated by background galaxies, (f) Highly significant 70 fim sources 
without 24 fim counterparts (S/N> 25). The gray outline indicates the approximate region covered by the combination of wavebands used 
to identify the different populations. 
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FlG. 12. — Color-color plots in the Spitzer and 2MASS wave- 
bands. The c olor of the symbols is identical to that employed in 
Figures [51 and [T0*1 Panels a . b, an d c, show plots similar to those 
modeled by W hitney et alJ 120031) for a variety of sources. The 
group 1 sources (green symbols) fall into the locii o f Class to 
Class II YSOs. Panel d, a plot similar to Figure 1 of Slo an et al.l 
( 2006 ) , clearly shows the two sequences of silicate and carbon rich 
dusty evolved stars (red symbols). 



colors are B-V n -0.20±0.12 and V-I » -0.25±0.07 
(after removing the line-of-sight reddening to the SMC) , 




Wavetenglh (urn) 

FlG. 13. — Average spectral energy distributions of the different 
populations described in the main text. The SEDs were obtained 
as the median at each wavelength over the entire distribution of 
objects. The drop in the SED at 4.5 (im seen in carbon stars and 
RSG and OAGB stars is probably due to the CO fundamental 
bandhead at 4.7 /an seen in absorption. 



while the photospheric colors of BO stars of metallicity 
Z = 0.1 - O.3Z are B — V = -0.27 and V - I = -0.36. 
The preferred geometry would be a face-on disk if the 
material is close to the star, a tenuous envelope, or a 
thin disk with a large central gap seen at arbitrary view- 
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ing angles. It could be argued that, since these stars 
were selected from their optical colors, the fact that they 
are only slightly reddened is not surprising. The FWHM 
of the 24 /im emitting B star sequence seen in Figure 
however, is only 0.08 mag in V — I. This is com- 
pletely consistent with the FWHM of the normal, main 
sequence B stars (0.07 mag). The 24 /mi-bright B stars 
therefore form a sequence as sharp as the main sequence 
with no sign of additional spread due to intrinsic redden- 
ing. Both this fact and Figure [7\ show that there is no 
evidence of a population of more reddened but otherwise 
similar sources. This suggests that the dusty material 
is arranged in a very tenuous envelope or in a thin disk 
with a large central gap, so that viewing angle does not 
play a role in the selection of these stars. Given the 
very large luminosity of the central source, material at 
Td«st ~ 150 K near an early-B st ar should be located 
at di stances of ~ 350 AU or more (Bac kman fe Parescd 
1993). 

What are these "24 /mi excess" early-B stars? Searches 
for these objects in the SIMBAD and VizieR databases 
indicate that out of approximately 190 stars in this 
group, 10 are eclipsing binari es (out of 1351 such object s 
found by OGLE in the SMC: iWvrzvkowski et al1l200l 
and 1 2 are emission line stars ijMevssonnier fc Azzopardl 
1993). In these aspects, this population does not seem 
peculiar. These sources are present in the SMC Bar and 
the tip of the Wing (FigurelTlTO , in regions of recent star 
formation. Their spatial distribution is, in fact, similar to 
that of the YSO candidates although not as strongly clus- 
tered, which is consistent with expectations for < 10 Myr 
old stars. We consider four possible explanations for the 
origin of these objects and their 24 /<m emission: 1) an 
accretion disk remnant, 2) free-free emission associated 
with emission line stars, 3) a circumstellar debris disk, 
or 4) nearby cirrus emission unrelated to the stars them- 
selves. We briefly discuss these possibilities below, but 
an in-depth discussion of these sources is deferred to a 
forthcoming paper (Bolatto et al., in prep.). 

Far infrared excess in stars is frequently a sign of stellar 
youth. These dusty B stars may be in the process of shed- 
ding their accretion remnants through photoevaporation. 
These early-B stars exhibit small ratios of IR excess lu- 
minosity to stellar luminosity, 10~ 4 < Lir/L* < 10~ 2 , 
rather than the much larger ratios observed in young 
stellar objects, Lir/Li* > 10 _1 . This fact, and their 
optical colors, suggest that their 24 /tm emission may 
be due to a remnant accretion disk or envelope with a 
large central gap, and that we may be observing them 
in the very late stages of the accretion process. Disk 
dispersal timescales, and in particular the timescale for 
removing the outer disk, are not well known for massive 
stars. Models of less massive stars, however, suggest that 
the outer disk clearing time is dominated by photoevap- 
oration, and that its dis persal occurs very quickl y after 
the inner disk is cleared l) Alexander et al.ll2006albT ). This 
phase in the lives of stars would be, consequently, very 
short, which suggests its observation is improbable. 

Although Milky Way emission line stars frequently 
have IR excesses, it seems unlikely that the 24 fim emis- 
sion observed in our SMC objects is explained by the 
various types of Be activity. A type of infrared excess 
associated with this activity is observed in Galactic Her- 
big Ae/Be stars, which are young, massive stars that 



may host circumstellar disks or envelopes that remain 
from the protostellar accretion ph ase and give rise to IR 
emission (Hill enbrand et al"]ll992j) . Classical Be stars, 
on the other hand, are older objects whose IR excess 
comes from free-free emission associated with their stel- 
lar winds and circumstellar disks. Be stars in the Milky 
Way tend to be later spectral types, but there is some ev- 
idence pointing to a large abundance of very early spec- 
tral types in the Magellanic C loud l)Orebel et al.lll993l 
iWisniewski fc Biorkmanl 12006(1 . We already mentioned 
that only a few of these 24 /mi emitting early-B stars are 
identified as Be stars in existing catalogs. Be envelopes 
are hot and clos e to the star, showing an i nfrared excess 
at 1-2 /mi fe.g.. lWaters fc Waelkensll998j) . Most, albeit 
not all, of our objects do have a small excess shortwards 
of 24 fim. At 4.5 fim the measured flux is 75% over 
the photosphere predicted by a Kurucz model tied to 
the observed V magnitude, and only a handful of these 
stars show an excess over a factor of 2. By compari- 
son, in the prototypical Herbig Ae/Be star AB Aur the 
observed emission exceeds the photospheric emission by 
more than an order of magnitude at 4.5 /im, and sim- 
ilarly large excesses are observed in m any classical Be 
stars (e.g., iMiroshnichenko et al .ll2003|) . Furthermore, 
the mean SED of our 24 /mi emitting B stars shortward 
of 24 fim is indistinguishable from the mean SED of B 
stars with the same V and V — I but undetected at 24 
/mi (Figure IT3l . More importantly, Figure fHfl (as well as 
inspection of the individual SEDs) shows that for most 
of our objects the observed short-wavelength excess can- 
not explain their 24 /im emission, where the SED takes 
a significant turn upwards. 

The low IR luminosity of these SMC stars in compar- 
ison to their photospheric stellar luminosit y is charac- 
teristic of systems hosting debris disks (e.g . jRieke et alJ 
l2005HGorlova et al.ll2006tlChen et aHl2006li . These disks 
are mostly devoid of gas, and form and replenish through 
the grinding of planetesimals. The small 24 /mi excess 
in our SMC early-B stars relative to their luminosity 
(10~ 4 < Lir/L* < 10~ 2 ) is very similar to the typical 
excess observed in debris disks (Lm/L* < 10 -3 ), where 
only a small fraction of the radiation from the central 
sourc e is captured by th e dust and reradiated in the IR 
(e.g.. lUzpen et ai1 l2005). This is, of course, an exciting 
possibility since it would provide indirect evidence for 
the existence of planets in a galaxy other than our own. 
There is a well-known deficit of giant planets (found by 
radial velocity techniques) around subsolar-metallicity 
stars that suggests such planets woul d be rare in a low- 
met a llicity galaxy such as the SMC l)Fischer fc Valentl 
2005). Recent observations, however, find that this 
deficit of giant planets does not translate into a defici t 
of debris disks (|Greaves et alJl2006t iBrvden et alj f2006) . 
Debris disks around massive early-B stars would be pop- 
ulated by millimeter-size particles, as radiation pressure 
should rapidly clear out smaller particles (unless gas drag 
is significant). By comparison, recent spectroscopic ob- 
servations suggest that debris disks around less massive 
A and F star s are populated by particles with siz es larger 
than 10 /mi ITura et al.H200l IChen et al.ll20f)l . 

The final possibility we consider for the origin of the 
24 /im emission is cirrus hotspots, sometim es noted 
as the Pleiades effect (e.g., iSloan et alJ 120041 . Cirrus 
hotspots occur where radiation from a massive star domi- 



19 



nates over the ISRF, heating a nearby patch of interstel- 
lar d ust above the ty pical tempe rature of cirrus emis- 
sion llvan Burenl I1989T) . Following Ivan Buren fc McCravl 
(1988), we estimate that the distance over which an 
early-B star can heat small (~ 0.1 /im size) ISM dust 
grains to temperatures Tdust ~ 150 K is r ~ 0.01—0.3 pc, 
or r < 1" at the distance of the SMC. This would ap- 
pear as a unresolved source in our data for everything 
except the brightest BO-09 stars, and thus be very dif- 
ficult to disentangle from a remnant accretion disk or 
a debris disk without additional observations at higher 
resolution. 

6. SUMMARY AND CONCLUSIONS 

We have used the Spitzer Space Telescope to image the 
star-forming body of the SMC in the mid and far IR with 
unprecedented sensitivity and resolution. 

IRAC images of the SMC show a wealth of diffuse emis- 
sion across the SMC seen for the first time. Particularly 
striking is the diffuse emission at 8 /im, often associated 
with star-forming regions. We have used these images 
to take a first look at the Fg/F24 ratio (an estimator of 
PAH abundance) across the SW Bar region of the SMC. 
We find that this ratio has large spatial variations, and 
takes a wide range of values even in a small region. Ex- 
tended emission in the southernmost part of the SW Bar 
has F 8 /F 2 4 > 0.7, typical of solar-metallicity galaxies, 
while near the sites of massive star formation and bright 
24 /im emission (suggesting large radiation fields) the ra- 
tio drop s to the Fs/F 2 4 < 0.1 o bserved in metal-poor 
galaxies (|Engelbracht et aTll2005t) . Because existing sur- 
veys of abundances in H II regions find no evidence for 
a gradient (or even a large scatter) for metallicity in the 
SMC, we conclude that the wide range of values taken 
by the Fs/F24 ratio in this galaxy, and the wide range 
of PAH abundance suggested by it, is primarily driven 
by the local ISRF and not metallicity. We find that the 
relationship between the Fg/F24 ratio and the FIR dust 
color represented by the F24/F70 ratio can be described 
by log[F 8 /F 24 ] « -0.71og[F 24 /F 7 o] - 1.4, showing that 
regions of blue FIR color are deficient in 8 /Ltm PAH emis- 
sion relative to their 24 /ini intensity. 

We combine these observations with optical and near- 
IR photometry in the literature to produce a catalog of 
point sources that spans from optical wavelengths to 70 
/im. We use this catalog to investigate the nature of the 
sources detected in the Spitzer wavebands. We find that 
the most prominent sources detected at the long wave- 
lengths fall into four categories: 1) sources with very faint 
optical counterparts and very red ([5.8] — [8.0] > 1.2) 
colors. This population corresponds to YSOs, with some 
contamination from background galaxies and foreground 
dwarfs that we estimate using the SWIRE ELAIS-N1 cat- 
alog. 2) Bright MIR sources (-12 < M 5 .$ < -7.5) with 
midly red colors (0.16 < [5.8] - [8.0] < 0.6), many of 
which have bright optical counterparts. These are iden- 
tified as carbon stars. 3) Bright mid infrared sources 
with neutral colors and bright optical counterparts, cor- 
responding to oxygen-rich AGB and red supergiant stars. 
And 4) a class of dusty early B stars (B3 to 09, based 
on the optical photometry) that are not appreciably red- 
dened optically, but that have large excesses at 24 /im 
most of which cannot be explained by free-free emission. 
We use our multiwavelength data to compile a catalog of 



282 bright YSOs in the SMC, where we choose the color- 
magnitude cuts to allow for only 7% contamination by 
background or foreground sources (Table 4) . 

Finally, we discuss the population of early B stars 
detected at 24 fan, which constitutes a few percent 
(~ 3 — 5%) of the B stars in a narrow region of the 
optical color-magnitude space, and shares many of the 
observational signatures of debris disks or cirrus hotspots 
in the Milky Way. Based on their small near IR excess 
and the fact that most of these objects are not found in 
existing objective prism catalogs, we conclude that the 
bulk of these systems are not emission line stars. They 
may constitute a population of young massive stellar sys- 
tems caught in the act of photoevaporating their accre- 
tion disks, but the timescale derived by models for such 
process suggests that this is unlikely. Dust associated 
with debris disks or cirrus hotspots (i.e., the Pleiades 
effect) appear to be the most likely explanation for the 
origin of their 24 /im emission. 

The S 3 MC point source catalog constitutes a deep and 
uniform photometric database for a large sample of stars 
located at the same distance, with similar metallicities, 
and low levels of foreground and background confusion. 
We expect that one of the contributions of this work 
will be to supply the basis on which to test different 
schemes for identifying mid- and far-infrared sources. 
The S 3 MC observations provide a unique resource for the 
study of stellar evolution, dust production, dust proper- 
ties, and the interaction between the ISM and star forma- 
tion in an environment that, because of its proximity and 
low metallicity, has closer resemblance to the conditions 
prevalent in primitive galaxies than any other galaxy that 
can be studied in similar detail. 
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